Search for: All records

The mixing of tracers by mesoscale eddies, parameterized in many ocean general circulation models (OGCMs) as a diffusive‐advective process, contributes significantly to the distribution of tracers in the ocean. In the ocean interior, diffusive contribution occurs mostly along the direction parallel to local neutral density surfaces. However, near the surface of the ocean, small‐scale turbulence and the presence of the boundary itself break this constraint and the mesoscale transport occurs mostly along a plane parallel to the ocean surface (horizontal). Although this process is easily represented in OGCMs with geopotential vertical coordinates, the representation is more challenging in OGCMs that use a general vertical coordinate, where surfaces can be tilted with respect to the horizontal. We propose a method for representing the diffusive horizontal mesoscale fluxes within the surface boundary layer of general vertical coordinate OGCMs. The method relies on regridding/remapping techniques to represent tracers in a geopotential grid. Horizontal fluxes are calculated on this grid and then remapped back to the native grid, where fluxes are applied. The algorithm is implemented in an ocean model and tested in idealized and realistic settings. Horizontal diffusion can account for up to 10% of the total northward heat transport in the Southern Ocean and Western boundary current regions of the Northern Hemisphere. It also reduces the vertical stratification of the upper ocean, which results in an overall deepening of the surface boundary layer depth. Finally, enabling horizontal diffusion leads to meaningful reductions in the near‐surface global bias of potential temperature and salinity.

We present idealized simulations to explore how the shape of eastern and western continental boundaries along the Atlantic Ocean influences the Atlantic meridional overturning circulation (AMOC). We use a state-of-the art ocean–sea ice model (MOM6 and SIS2) with idealized, zonally symmetric surface forcing and a range of idealized continental configurations with a large, Pacific-like basin and a small, Atlantic-like basin. We perform simulations with five coastline geometries along the Atlantic-like basin that range from coastlines that are straight to coastlines that are shaped like the coasts of the American and African continents. Changing the Atlantic basin coastline shape influences AMOC strength in a manner distinct from simply increasing basin width: widening the basin while maintaining straight coastlines leads to a 10-Sv (1 Sv ≡ 10⁶m³s⁻¹) increase in AMOC strength, whereas widening the basin with the geometry of the American and African continents leads to a 6-Sv increase in AMOC strength, despite both cases representing the same average basin-width increase relative to a control case. The structure of AMOC changes are different between these two cases as well: a more realistic basin geometry results in a shoaled AMOC while widening the basin with straight boundaries deepens AMOC. We test the influence of the shape of the both boundaries independently and find that AMOC is more sensitive to the American coastline while the African coastline impacts the abyssal circulation. We also find that AMOC strength and depth scales well with basin-scale meridional density difference, even with different Atlantic basin geometries, illuminating a robust physical link between AMOC and the North Atlantic western boundary density gradient.

We document the configuration and emergent simulation features from the Geophysical Fluid Dynamics Laboratory (GFDL) OM4.0 ocean/sea ice model. OM4 serves as the ocean/sea ice component for the GFDL climate and Earth system models. It is also used for climate science research and is contributing to the Coupled Model Intercomparison Project version 6 Ocean Model Intercomparison Project. The ocean component of OM4 uses version 6 of the Modular Ocean Model and the sea ice component uses version 2 of the Sea Ice Simulator, which have identical horizontal grid layouts (Arakawa C‐grid). We follow the Coordinated Ocean‐sea ice Reference Experiments protocol to assess simulation quality across a broad suite of climate‐relevant features. We present results from two versions differing by horizontal grid spacing and physical parameterizations: OM4p5 has nominal 0.5° spacing and includes mesoscale eddy parameterizations and OM4p25 has nominal 0.25° spacing with no mesoscale eddy parameterization. Modular Ocean Model version 6 makes use of a vertical Lagrangian‐remap algorithm that enables general vertical coordinates. We show that use of a hybrid depth‐isopycnal coordinate reduces the middepth ocean warming drift commonly found in purez^*vertical coordinate ocean models. To test the need for the mesoscale eddy parameterization used in OM4p5, we examine the results from a simulation that removes the eddy parameterization. The water mass structure and model drift are physically degraded relative to OM4p5, thus supporting the key role for a mesoscale closure at this resolution.